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SUMMARY 

A  fixed  based  simulator  study  of  the  approach  and 
tangential  landing  of  a  class  of  manned  re-entry  gliders 
was  conducted  to  determine  the  piloting  techniques  and  ve- 
hicle parameters  which  best  define  the  approach  and  land- 
ing problem.  Tests  were  conducted  on. 5  configurations 
typical  of  proposed  vehicles  with  maximum  lift  to  drag 
ratios  of  3  to  5»   The  problem  was  Investigated  from  final 
approach  to  touchdown,  In  still  and  turbulent  air,  simulat- 
ing longitudinal  equations  of  motion  only. 

A  three  phase  landing  consisting  of  a  steady  glide 
final  approach,  a  nearly  constant  load  factor  flare  maneuver, 
and  a  15  to  20  second  shallow  angle,  final  glide  was  the 
preferred  technique.   The  vehicle  parameters  found  to  be 
most  Important  were  maximum  lift  to  drag  ratio  and  speed 
for  minimum  drag.   A  landing  analysis  based  on  fitted  para- 
bolic* drftg  pol«*rs  A«4  ft  ponpUnt  3-oaO.  faofcor  T\iM*<*    was 
found  to  be  valid  for  the  vehicles  tested. 


INTRODUCTION 

The  manned  space  vehicle  capable  of  atmospheric  re- 
entry and  aerodynamic  maneuvering  to  a  horizontal  landing 
is  presently  conceived  as  a  fixed  geometry  wing-like  lift- 
ing body  of  low  aspect  ratio  and  high  profile  drag.   The 
vehicle  may  have  relatively  low  wing  loading  and  no  pro- 
pulsion in  the  landing  phase. 

The  approach  and  landing  of  such  vehicles  is  character- 
ized by  high  approach  speeds,  steep  approach  angles,  rapid 
deceleration  in  near  horizontal  flight,  and  the  probability 
of  high  speed  and  sink  rate  at  touchdown.   These  character- 
istics pose  a  piloting  task  which  demands  concentration 
and  accuracy,  while  leaving  scant  margin  for  error.   This 
is  particularly  true  if  landing  facilities  and  navigation- 
guidance  aids  similar  to  those  which  accommodate  modern 
high  performance  fighter  aircraft  are  to  be  used. 

Approach  and  landing  techniques  which  minimize   the 
pilot  judgment  required,  avoiding  the  usual  flare  or  pull- 
up  just  prior  to  touchdown  and  providing  reasonable  capa- 
bility of  pilot  escape  with  today's  ejection  equipment, 
are  proposed  in  REFS.  1,  2,  3,  and  5»   The  technique  has 
been  used  successfully  with  the  X-15  research  airplane 
(REF.  k). 


The  objective  of  this  research  was  to  study  the  tan- 
gential landing  of  a  class  of  re-entry  gliders  using  a 
fixed-base  simulator  with  emphasis  on  determining  the 
pilot  and  vehicle  parameters  which  best  define  the  approach 
and  landing  problem. 

Simulated  landings  of  vehicles  with  maximum  lift  to 
drag  ratios  of  3  to  5  were  conducted  in  still  and  turbulent 
air.   The  landing  performance  was  analyzed,  and  an  analytical 
procedure  applicable  to  generalized  vehicles  developed. 

The  principle  pilot  and  author  was  a  designated  Naval 
Aviator  with  eight  years  of  military  flight  experience  and 
over  2^00  hours  In  propeller  type  airplanes. 


EQUIPMENT  AND  PROCEDURE 

1.   Simulator. 

The  simulator  consisted  basically  of  a  fixed  base  seat, 
Instrument  panel,  CRT  runway  display,  and  side  controller 
(FIG.  1).   The  computing  equipment  Included  an  Electronic 
Associates  Pace  analog  computer,  a  Llbratrol  50°  digital 
computer,  with  an  analog-digital  voltage  converter.   Data 
recording  equipment  Included  a  six  channel  Vlsicorder  906C 
strip  recorder  and  the  analog  X-Y  plotter.   A  detailed 
description  of  the  simulator  and  associated  equipment  is 
contained  in  Appendix  A. 

The  pilot  was  provided  instrument  display  of  the  ve- 
hicle's altitude,  airspeed,  pitch  attitude,  flight  path 
angle,  and  glide  path  error.   The  runway  simulation  con- 
sisted of  a  porthole  view  of  the  runway  boundaries  shown 
in  standard  perspective.   The  simulated  dimensions  of  the 
runway  were  10,000  feet  long  by  200  feet  wide.   Airplane 
pitch  altitude  cues  were  represented  by  vertical  motion 
and  position  of  the  runway  display.   Visual  range  and 
height  cues  could  be  interpreted  from  the  perspective 
view  of  the  landing  area.   Only  motion  in  the  plane  of 
symmetry  was  simulated  and  no  lateral-directional  task 
was  assigned  the  pilot. 


The  side  controller  was  given  essentially  constant 
force  versus  displacement  characteristics  with  no  centering 
forces-   Breakout  forces  and  dead  zone  were  negligible. 
No  control  system  dynamics  were  simulated. 

Atmospheric  turbulence  was  simulated  by  introducing  a 

rectangular  or  uniform  probability  density  of  angle  of 

attack  (<^)  disturbances  into  the  airplane  equations.   The 

disturbances  could  be  converted  to  an  equivalent  vertical 

gust  velocity  by      ^  _  y, 

Y  7 

A  vertical  gust  of  6.5  fps  RMS  at  the  speed  for  minimum 

drag  was  used. 

2.   Vehicle  Dynamics. 

The  basic  longitudinal  equations  of  motion  for  a  rigid 
body  with  the  standard  axis  system  fixed  in  the  aircraft 
are  (REF.  6): 

(1)  jl-  Wa^Of-  ™VZt  ~°  (Lift) 

(2)  i)-/M^o^/-  ^  Jt  ^  °  (Drag) 

( 3 )  /Sy  -  Ty&  -■  o  ( Moment ) 
The   kinematic  equations  are: 

(Li.)  ^  ^   i/s<u~^  cK  =    l/c<rf<K  A ^  ■/  i^u~t  rQ 

(6)  <jj  -   4  -  (*-y<p )  ~fa*^> ^/ 


where  the  terms  sin  ar  and  cos  fr     have  been  linearized 
using  a  perturbation  about  X,    ,  a  nominal  flight  path  angle. 

The  aerodynamic  forces 

/r-  /  («,  *) 

were  expanded  In  a  Taylor  Series  Including  second  order 
terms  where  necessary  to  represent  the  re-entry  vehicle 
aerodynamics. 

(7)    t- &*/*.+&+  3*1 

(9)  /7-jfrJAi.JC^  +    ^   «  +  2z    *  +  2«&         \ 

Aerodynamic  forces  were  assumed  Independent  of  Mach 
number  for  this  experiment  since  the  maximum  Mach  number 
simulated  was  approximately  .50.   Experiments  with  higher 
Mach  numbers  should  consider  compressibility  effects. 

Equations  1  through  6  with  7,  8,  and  9  were  mechanized 
to  simulate  the  vehicle  dynamics  in  the  experiment.   A 
detailed  description  of  the  programming  and  the  computa- 
tion of  the  runway  display  are  given  in  Appendix  A. 


3.  Vehicle  Characteristics. 

In  order  to  evaluate  the  possible  effects  of  re-entry- 
vehicle  non-llnearltles,  It  was  decided  to  fit  as  accurately 
as  possible  the  wind  tunnel  data  of  typical  vehicles,  B  and 
C  (FIGS.  2  and  3).      Tests  were  also  conducted  on  a  slightly 
simplified  vehicle  designated  B,  In  which  the  second  order 
aerodynamic  coefficients  were  eliminated.   Two  additional 
hypothetical  vehicles  (A  and  D)  of  the  same  general  class 
were  tested.   In  each  test  the  vehicle  was  programmed  to 
simulate  the  desired  characteristics  in  both  a  range  of 
trimmed  and  untrlmmed  flight  conditions  most  appropriate 
to  the  individual  glider. 

The  programmed  characteristics  of  the  five  vehicles 
tested,  A,  B,  B  ,  C,  and  D,  are  shown  in  TABLE  I  and  FIG.  ^. 
The  characteristics  for  all  vehicles  tested  are  in  the 
clean  configuration.   The  lowering  of  landing  gear  and 
flaps  was  not  simulated  In  this  investigation. 

4.  Simulator  Flight  Task. 

The  flight  task  was  designed  to  simulate  a  straight 
in  approach  and  landing  on  a  10,000  foot  runway  using 
essentially  standard  aircraft  instruments  and  visual  dis- 
play.  In  preliminary  experiments  the  problem  was  started 
with  the  vehicle  at  20,000  feet  altitude.   However,  the 
approach  phase  with  the  aircraft  Initially  positioned  on 
a  nominal  approach  path  was  found  to  be  trivial  and  the 
test  runs  were  shortened  to  save  computer  time. 


Each  test  run  was  commenced  In  an  equilibrium  glide 
condition,  initially  positioned  on  a  nominal  glide  path 
to  the  runway  at  an  altitude  of  5,000  feet.   The  initial 
velocity  and  flight  path  angle  were  varied  over  the  range 
of  practical  flight  conditions. 

In  each  approach  configuration,  the  altitude  at  which 
the  flare  or  pull-up  was  initiated  was  varied  until  the 
landing  performance  for  that  configuration  was  considered 
optimum  by  the  pilot.   A  series  of  at  least  5  landings 
for  each  approach  condition  was  then  made  in  which  the 
pilot  attempted  to  execute  the  best  possible  approach, 
flare,  and  landing  maneuver. 

The  glider  trajectory  on  all  test  runs  was  recorded 
by  plotting  V  versus  Y  on  the  X-Y  plotter.  Five  of  the 
following  variables  were  normally  recorded  as  a  function 
of  time: 

(a)     V*  {J)     A 

tij    *  a)   <£ 

Cc)      X        (f)      0 


Pilot  comments  were  made  and  recorded  on  each  test  run. 
No  minimum  speed  for  touchdown  was  specified.   The 
maximum  speed  for  touchdown  was  specified  by  an  assumed 
maximum  vertical  speed  of  10  fps.   Landing  performance 


8 


was  Judged  on  controllability  and  accuracy  for  the  flare 
maneuver,  touchdown  accuracy,  ability  to  maintain  vertical 
speed  at  touchdown  within  assumed  limits,  and  overall 
consistency  of  the  landing  maneuver.   Additional  comments 
were  made  on  the  pilot's  opinions  of  the  short  period  ve- 
hicle dynamic  characteristics  where  applicable. 

Simulated  turbulence  with  an  RMS  equivalent  to  6.5  fps 
for  the  vehicle  average  velocity  during  final  glide  was 
used  In  tests  of  vehicles  A,  B,  and  D.  Pilot  comments  on 
the  effects  of  turbulence  were  also  made. 


RESULTS  AND  DISCUSSION 

It   Description  of  the  Landing  Maneuver. 

An  evaluation  of  the  preliminary  experiments  validated 
the  description  .(REFS.  1  through  5)  of  the  landing  maneuver 
as  divided  into  five  parts  or  phases. 

The  first  part  is  a  pilot  guidance  phase  probably 
beginning  after  re-entry  during  which  the  pilot  maneuvers 
to  a  high  key  position  near  or  over  the  landing  field  aided 
by  a  radar-voice  link  and/or  Inert ial  navigation  systems. 

The  second  part  is  usually  described  as  a  visual 
straight  in  or  circling  pattern  in  which  the  glider  is 
maneuvered  to  intercept  the  final  approach  leg  at  a  pre- 
scribed altitude  and  airspeed.   During  this  maneuver  the 
pilot  might  use  TACAN  equipment  as  a  primary  navald. 

The  third  and  subsequent  parts  are  the  subject  of 
this  report  and  will  hereafter  be  referred  to  as  PHASES  I, 
II,  and  III  (FIG.  5). 

The  three  phases  are  described  as  follows: 

PHASE  I  -  FINAL  APPROACH.   The  final  approach  is  a 
steady  state  glide  at  nearly  constant  airspeed  and  flight 
path  angle  aimed  at  a  fixed,  predetermined  geographic  point 
short  of  the  runway. 

PHASE  II  -  FLARE.   The  flare  is  essentially  a  constant 
load  factor  pull-up  initiated  at  a  predetermined  altitude 
during  which  the  sink  rate  is  reduced  to  near  the  maximum 
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allowable  for  touchdown.   The  maneuver  Is  completed  at  a 
typically  low  altitude  and  small  flight  path  angle. 

PHASE  III  -  FINAL  GLIDE.   The  aircraft  is  flown  in 
essentially  one  "g",  constant  ¥   flight  as  excess  speed 
and  altitude  are  reduced  for  touchdown.   During  this  phase 
the  pilot  has  rather  precise  control  over  his  altitude  and 
flight  path  angle  but  essentially  no  control  over  decelera- 
tion.  Landing  gear  and  flaps  are  normally  lowered  in  this 
phase. 

Another  phase  which  Is  not  included  in  this  work  is 
the  aircraft  roll-out  or  slide-out  after  runway  contact. 

2.   Quality  of  Simulation. 

The  preliminary  evaluation  of  the  simulator  Indicated 
that  the  combination  of  instrumentation  and  visual  display 
provided  sufficient  realism  to  enable  a  pilot  to  perform 
and  evaluate  the  final  approach  and  landing  maneuver.   The 
absence  of  depth  perception  in  the  runway  display  was  an 
early  difficulty  solved  .by  locating  the  flight  path  angle 
Indicator  adjacent  to  the  CRT  display  where  it  could  be 
Included  as  a  peripheral  supplement  to  the  limited  physical 
depth  cues.   Accurate  landings  could  then  be  performed  with 
reference  to  the  scope  for  attitude  cues,  and  with  airspeed, 
altimeter,  and  flight  path  indicator  providing  precise  con- 
trol of  sink  rate. 
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During  the  flare,  the  flight  path  indicator  and  altimete: 
provided  the  primary  cues  for  controlling  the  flight  path 
curvature  in  such  a  way  as  to  arrive  at  the  desired  h  and  X* 
for  final  glide.   The  CRT  display  and  attitude  indicator 
provided  the  necessary  attitude  cues. 

3.   Selection  of  Touchdown  Speed. 

Since  the  final  approach  and  landing  problem  can  be 
divided  into  3  rather  distinct  phases,  the  analysis  may 
be  divided  in  the  same  manner  by  simply  matching  boundary 
conditions.   It  is  also  convenient  to  start  the  problem 
in  reverse  since  it  was  found  that  the  touchdown  conditions 
can  be  determined  independent  of  the  approach,  flare,  and 
final  glide. 

In  general  the  selection  of  a  minimum  ground  contact 
speed  is  not  dependent  on  a  conventional  "stall  speed". 
The  re-entry  vehicle,  because  it  is  most  probably  a  lift- 
ing body,  has  no  well  defined  Cr  maximum.   Instead  the 
minimum  touchdown  speed  might  be  determined  by  any  one  of 
the  following: 

a.  Wing  buffet 

b.  Directional  stability  limits 

c.  Ground  clearance  limits 

d.  Maximum  landing  gear  extension 

e.  Minimum  trim  speed 
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For  this  experiment  It  was  assumed  that  there  were 
no  limitations  on  the  minimum  speed  except  the  ability  of 
the  pilot  to  precisely  control  the  touchdown.   Elevator 
control  power  was  provided  in  excess  of  that  required  to 
trim  at  the  touchdown  speed.   Pitch  damping  was  Incorpo- 
rated in  all  vehicles  to  insure  satisfactory  short  period 
dynamic  characteristics  at  minimum  speed.   Since  the  long 
period  vehicle  dynamics  are  of  paramount  interest  in  this 
study,  difficulties  with  short  period  motions  were  thus 
virtually  eliminated,  (See  Section  10  on  Short  Period 
Characteristics).  Under  these  conditions  the  vehicle 
could  be  flown  to  a  touchdown  at  any  speed  within  the 
acceptable  limits  of  vertical  speed  at  touchdown  assumed 
(10  fps). 

It  was  found  in  the  experiment  that  if  flight  prog- 
ressed to  a  speed  below  a  value  characteristic  of  each 
vehicle,  the  loss  of  kinetic  energy  was  so  rapid  that  an 
excessive  sink  rate  usually  resulted.   This  indicated  that 
the  pilot  was  unable  to  increase  the  angle  of  attack  rapidly 
enough  to  maintain  one  Mg"  flight  (L=W).   This  characteris- 
tic speed  was  found  to  be  slightly  below  the  speed  for 
minimum  drag  (Vg). 

Assuming  a  characteristic  parabolic  drag  polar  where 
(10) 


^-  3b  *  ^ 


13 


the  CL  for  maximum  lift  to  drag  ratio  is 


Then 


4, 


*--#&£ 


and 

(13)  £■-£-   — 


Dm**."  ZiKCc^ 


Letting 


(15)  n-  £ 


Substituting   In 

(16) 

0-jfKt<\.    +fy> 

yields 
(17) 

*-jfr*"l 

Then  with  the  drag  equation  (2) 

(18)  <& £     <  jt±±*£*~r  +  M 
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2 
In  this  equation  the  variable  u  represents  a  linear  con- 
tribution of  the  profile  drag  to  the  deceleration  while 
1/u   is  a  hyperbolic  contribution  of  the  induced  drag. 
The  constant  terms  are  a  combination  of  gravitational 
acceleration  and  vehicle  characteristics. 

For  flight  speeds  less  than  Vg  the  induced  drag  pre- 
dominates and  deceleration  rapidly  approaches  infinity 
as  shown  in  FIG.  6. 

A  simple  statistical  analysis  of  the  successful 
simulator  landings  are  shown  in  FIG.  ?.   Based  on  this 
empirical  data  a  minimum  touchdown  speed  for  a  reasonable 
probability  of  success  should  be  approximately  .88  V^  in 
the  absence  of  other  requirements  for  a  higher  landing 
speed  as  previously  discussed. 

4.   The  Maximum  Touchdown  Speed. 

The  maximum  allowable  touchdown  speed  is  generally 
a  structural  limitation  imposed  by  the  vehicle.   The  ver- 
tical speed  at  ground  contact  depends  on  both  V  and  o 
since  <JA 

FIG.  8  shows  h  versus  &    where  the  limiting  sink  rate  is 
assumed  to  be  10  fps.   In  practice  the  airplane  will  be 
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landed  with  some  small  but  finite  2f   found  to  be  minus  2 
to  2.5  degrees  In  this  experiment.   FIG.  8  shows  a  dis- 
advantage of  vehicles  with  high  VR,  hence  high  landing 
speeds,  requiring  relatively  smaller  &*   at  touchdown. 
It  was  found  with  all  the  vehicles  flown  In  this 
experiment  that  a  maximum  sink  rate  of  10  fps  is  a  prac- 
tical limitation  and  imposes  no  piloting  problems  unless 
landings  are  attempted  at  speeds  below  .88  Vq.   These 
landings  usually  occurred  unintentionally  If  the  pilot 
leveled  the  vehicle  high  and  subsequently  ran  out  of  kinetic 
energy  before  he  could  lower  the  vehicle  to  the  ground. 

5.   Final  Glide. 

The  final  glide  (PHASE  III)  was  found  best  character- 
ized by  the  time  T  between  the  completion  of  PHASE  II  and 
touchdown.   This  "time-to-go"  to  touchdown  is  the  interval 
needed  by  the  pilot  to  correct  the  errors  in  Judgment  and/or 
execution  which  occur  during  flare  allowing  him  to  pre- 
cisely lower  the  rapidly  decelerating  glider  to  the  runway 
with  an  acceptable  sink  rate. 

Experience  with  the  X-15  airplane  (REF.  4)  and  flight 
simulator  experiments  with  fighter  aircraft  in  high  drag 
configurations  (REFS.  1,  2,  and  3)1  Indicates  that  a  min- 
imum of  15  to  20  seconds  should  be  allowed  In  the  final 
glide.   Reference  3  states  that  30  seconds  is  usually  the 
maximum  required  for  any  condition. 
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In  the  simulator  flights,  10  seconds  or  less  was  found 
generally  unsatisfactory  resulting  In  unacceptable  "scatter" 
In  landing  performance  measured  In  terms  of  touchdown  posi- 
tion, sink  rate,  and  speed  dispersion,   10  to  15  seconds 
was  generally  found  to  be  marginal  resulting  In  some  scatter 
particularly  In  turbulent  air.   15  to  20  seconds  was  found 
to  be  quite  satisfactory  for  all  vehicles.   Times  greater 
than  20  seconds  were  seldom  experienced  and  cannot  be 
evaluated. 

Since  the  minimum  touchdown  speed  is  determined  inde- 
pendent of  the  trajectory,  integration  of  the  drag  equation 
over  the  final  glide  allows  the  prediction  of  the  speed 
required  at  flare  completion  for  any  desired  deceleration 
time. 

Equation  (18)  In  the  form 
(19)  tfk  _    f  U*-t  J£*U~<  y  -h  /? 

where  <J  c/ 

expresses  the  final  glide  deceleration  for  parabolic  drag 
assumptions* 

Equation  (19)  integrated  Is 

a->-±  —        I/*?  r-e 


1? 


where 


G- i£*+i2*u  +  {f%/)~  fur*. 

a  +  fzu. 


Z4  /Z~i 


U- 


Computation  with  equation  (20)  Is  unhandy  and  a  simple 
graphical  solution  might  be  constructed  for  a  particular 
value  of  E  by  plotting  -jr.    backwards  from  touchdown  using 
a  nominal  <#  of  about  minus  2  degrees.   The  function  Is 

7 

quite  linear  In  the  range  of  interest  and  sufficient  accu- 
racy can  be  obtained  by  plotting  only  a  few  points.   FIG.  9 
is  the  approximate  solution  of  equation  (20)  for  the  ve- 
hicles tested  assuming  touchdown  at  .88  VR.   Results  of 
the  experiment  show  very  close  agreement  with  this  analysis 
tending  to  show  the  validity  of  the  parabolic  drag  polar 
assumption.   FIG.  10  shows  a  comparison  of  trimmed  aero- 
dynamic characteristics  for  vehicles  B  and  C  with  their 
parabolic  approximations.   For  these  examples,  the  region 
of  interest  for  flare  and  landing  lies  between  C^  of 
approximately  .15  and  .^0.   This  region  shows  quite  rea- 
sonable agreement. 
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The  deceleration  characteristics  of  re-entry  vehicles 
and  their  landing  speeds  define  the  problem  in  the  final 
glide  where  a  rather  fixed  time  must  be  allowed.   Further- 
more, the  deceleration  is  approximately  constant.   Consider 
then  the  form  of  the  drag  equation  for  final  glide 


Substituting  the  velocity  ratio  at  VR 


(22)  ^~_ 


Neglecting  the  gravitational  force  for  simplicity 
(23)  </y .  a 

72~  --* 


Thus  it  can  be  seen  that  the  vehicles  with  higher  lift  to 
drag  ratios  will  have  superior  deceleration  characteristics. 
Next  consider  (23)  Integrated  between  flare  completion  and 
touchdown. 


(24) 


7~ 


fc-  fe  +  "* 


rr>/fi/. 


But  if  Vmin  equals  .88  VR 

(25)  ^ _ 
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This  simplified  equation  shows  the  effect  of  vehicle  para- 
meters on  final  glide  and  landing  speed.  The  final  glide 
characteristics  are  necessarily  determined  by  both  decel- 
eration, a  function  of  E.  alone,  and  the  speed  for  minimum 
drag,  a  function  of  both  wing  loading  and  the  drag  char- 
acteristics.  An  example  of  this  is  shown  in  the  table 
below.   Values  shown  were  calculated  for  T  of  15  seconds 
using  FIG.  9. 


VEHICLE 

E 

w 
s 

K 

0 

Vr 

Vmin 

Vfg 

A 

3.0 

22.7 

.463 

.06 

230 

202 

317 

B 

5.0 

29.4 

.318 

.0316 

280 

246 

286 

C 

4.4 

25.0 

.381 

.034 

265 

233 

288 

D 

4.4 

31.4 

.258 

.050 

265 

233 

288 

Note:   K  and  CD  are  based  on  the    idealized  drag  polar. 

The  above  table  depicts  four  of  the  five  vehicles  tested. 
Vehicle  B-,  is  identical  to  B  except  in  the  computer  program- 
ming.  Vehicles  B  and  C  are  representative  vehicles  of  the 

are 
re-entry  class  while  A  and  D/ hypothetical  vehicles. 

It  may  be  noted  that  vehicles  C  and  D,  although  they 
have  different  vehicle  parameters,  show  identical  landing 
characteristics.   This  is  possible  since  VR  is  defined  by 
3  vehicle  parameters  while  E  is  defined  by  2.   Comparing 
vehicles  B  and  C  shows  the  effect  of  wing  loading.   Vehicle 
B  is  superior  in  maximum  lift  to  drag  ratio  but  must  land 
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at  the  higher  speed.   Comparing  B  to  A  shows  the  same 
tendency  but  even  more  graphically.   However,  vehicle  A 
must  be  provided  the  greater  starting  speed  and  distance 
from  landing  area.   From  the  simulator  study  It  could  not 
be  determined  what  the  effect  of  the  higher  deceleration 
and  greater  landing  speeds  would  be  on  pilot  opinion  of 
the  vehicle.   The  measure  of  success  seemed  to  depend 
primarily  on  the  T.   If  the  vehicle  could  complete  the 
flare  with  sufficient  kinetic  energy  for  approximately 
15  seconds  or  more  In  final  glide,  the  landing  could  be 
successfully  completed  independent  of  the  relative  decel- 
eration and  speed  characteristics. 

6.   Flare. 

The  flare  maneuver  (PHASE  II)  was  found  to  be  best 
characterized  by  the  normal  load  factor  used  during  the 
maneuver.   Although  it  was  found  that  the  flare  was  essen- 
tially a  "closed  loop"  maneuver  In  which  the  pilot  senses 
his  airplane's  flight  path  curvature  and  altitude,  the 
maneuver  tends  to  be  accomplished  at  nearly  constant  load 
factor.   In  this  experiment,  the  pilot  was  not  provided 
an  accelerometer,  yet  repeatedly  the  flare  found  best  liked 
in  terms  of  controllability  and  accuracy  was  about  1.2  to 
1*3  "g".   A  series  of  typical  time  histories  of  ^/£ 
calibrated  in  "g"  units  are  shown  In  FIG.  11. 
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Representing  "favorable"  or  near  "optimum"  pilot 
comments,  these  data  show  the  consistent  average  of  the 
1.2  to  1.3  "g"  flare  for  k   out  of  the  5  vehicles  over  a 
variety  of  Initial  velocity  and  altitude  conditions. 
Moreover,  the  simulator  tests  showed  that  the  favorable 
or  most  preferred  flare  load  factor  was  actually  Independ- 
ent of  altitude  and  velocity  of  Initiation,  the  pilot  tend- 
ing to  adjust  the  flare  Initiation  altitude  to  give  the 
preferred  load  factor  over  a  number  of  trial  runs,  and 
without  reference  to  an  accelerometer  or  the  data  of  that 
particular  series  of  test  runs.   Flares  using  higher  "g" 
tended  to  be  difficult  to  control  resulting  most  often 
In  misjudged  level-off  altitudes  and  subsequent  poor  land- 
ing or  even  catastrophic  over  shooting  and  premature  ground 
contact.   No  time  history  data  is  available  on  vehicle  A, 
but  the  same  tendency  was  evident  In  simulated  flight 
experience. 

Analysis  of  the  flare  maneuver  began  by  considering 
equation  (18) 

(26)        ^  _-£  f ^ *+ t^*"** *+«. y 

together  with  the  lift  equation  (1)  in  the  form 
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For  -»  -  M. 


equation  (2?)  is 

(28)  ^K.  fn-t) 

Then 

(29)  <&        /  *£("-')  u.  ? 


To  show  the  effect  of  the  constant  load  factor  assump- 
tion, a  series  of  trajectories  were  calculated  for  comparison 
with  actual  test  trajectories  in  V  versus  X*   space  as  re- 
corded on  the  X-Y  plotter  (FIG.  12).   It  may  be  noted  that 
the  n  equal  to  1.2  step  function  approximation  to  the  actual 
"closed  loop"  load  factor  profile  gives  excellent  results 
in  predicting  the  actual  trajectory. 

It  would  appear  that  the  higher  Hg"  trajectory  would 
be  superior  since  the  conservation  of  kinetic  energy  for 
use  in  the  final  glide  is  paramount.   However,  the  problem 
of  precisely  controlling  the  flight  path  curvature  at  the 
higher  turning  rates,  lends  favor  to  the  less  severe  flare. 

Equation  (29)  reveals  that  the  slope  of  the  trajectory 
in  V  versus  <K  space, which  is  related  to  the  very  important 
kinetic  energy  loss  in  the  flare,  is  a  strong  function  of  E. 
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To  Illustrate  this,  a  family  of  trajectories  for  vehicles 
A,  B,  and  C  were  constructed  using  the  load  factor  1.2  a 
consistent  with  the  empirical  data  (FIG.  13).   All  vehicles 
are  started  with  the  steady  glide  at  450  fps.   The  final 
glide  condition  In  each  case  Is  minus  2  degrees.   Equation 
(29)  was  used  and  the  non-dimensional  u  converted  to  the 
more  familiar  co-ordinate  V. 

It  may  be  observed  that  the  kinetic  energy  lost  during 
the  flare  increases  rapidly  as  E  decreases.   Since  the 
steady  glide  speed  Is  limited  by  compressibility  to  about 
kys   fps,  and  the  flare  completion  speed  Is  limited  by  the 
necessity  of  a  finite  deceleration  time,  there  Is  a  defi- 
nite limit  to  the  capability  of  a  vehicle  to  perform  the 
landing  maneuver  using  this  technique.   This  was  found  In 
test  experience  with  vehicle  A  wherein  the  T  final  glide 
was  found  to  approach  zero  and  landing  performance  was 
unsatisfactory.   Once  in  the  final  glide  the  lower  E  char- 
acteristics (which  might  be  typical  of  a  vehicle  with  gear 
and  flaps  extended)  does  not  seem  to  be  a  serious  limitation 
If  sufficient  initial  velocity  can  be  provided.   Use  of  a 
higher  load  factor  in  the  flare  decreases  the  kinetic  energy 
loss,  but  again  the  precision  required  seems  to  limit  this 
technique. 

Another  serious  vehicle  limitation  might  be  the  alti- 
tude lost  during  the  flare.   If  the  flare  must  be  Initiated 
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at  a  relatively  high  altitude,  the  accuracy  with  which  the 
pilot  can  perform  the  maneuver  must  suffer.   In  this  re- 
spect the  simulator,  even  with  Its  "Ideal"  altimeter  and 
limited  physical  cues,  showed  this  to  be  the  case.   As  the 
flare  altitude  Increased,  the  pull-up  maneuver  tended  to 
become  "open-loop"  until  much  of  the  altitude  was  dissipated 
At  some  lower  altitude  a  second  or  modified  flare  was 
usually  evident.   In  general, the  high  altitude  maneuver 
was  more  difficult  and  less  precise. 

To  calculate  the  altitude  loss  In  a  flare  consider 
(30) 


or 
(3D 


K*usrc<>r 

With  equation  (28) 
(32) 


This  calculation  may  be  carried  out  stepwise  together  with 
the  trajectory  calculation.   However,  by  assuming 


^£~   4*\i-  ccA/sr. 
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Then 


03)  7     ._*    3 


3<?(n-t) 


*~t-  av 


With  equation  (33)  the  altitude  lost  In  the  flare  for  the 
three  trajectories  of  FIG.  12  is 

VEHICLE  A  4780  ft. 
VEHICLE  B  1398  ft. 
VEHICLE  C       1620  ft. 

These  results  are  consistent  with  the  experiment  and  accurate 
enough  for  planning  an  approach  since  the  pilot  will  actually 
vary  his  pull-up  considerably  depending  on  his  Judgment. 

Simulator  experience  with  the  flare  maneuver  through 
more  than  about  2,000  ft.  showed  them  to  be  difficult  to 
perform  precisely.   Flare  maneuvers  at  less  than  about  600  ft. 
also  proved  to  be  undesirable  but  these  are  seldom  encountered 
if  the  1.2  Mgw  pull-up  is  used. 

7.   Final  Approach. 

The  steady  glide  flight  conditions  of  the  final  approach 
phase  (PHASE  I)  can  be  determined  by  calculating  a  desired 
trajectory  backwards  from  touchdown  until  an  equilibrium 
glide  flight  condition  is  reached  for  any  given  vehicle. 

The  speed,  flight  path  angle,  and  rate  of  descent  In 
the  final  approach  will,  of  course,  depend  on  the  lift-drag 
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characteristics  of  the  vehicle.   However,  the  speed  Is 
bounded  by  compressibility  effects  and  the  maximum  accept- 
able rate  of  descent  probably  depends  on  pilot  opinion. 
The  latter  limitation  was  not  a  difficulty  In  the  simulator 
and  therefore  could  not  be  evaluated* 

8.   Synthesis  of  an  Approach. 

An  example  of  how  an  approach  could  be  synthesized 
Is  given  here  and  shown  In  FIG.  lb.      Vehicle  C  Is  used  for 
this  example. 

It  was  assumed  that  the  touchdown  speed  was  not  other- 
wise limited  and  .88  VR  was  used.   It  was  further  assumed 
that  15  seconds  was  desired  In  final  glide,  and  the  lift- 
drag  characteristics  were  constant  (no  gear  and  flaps 
lowered  In  PHASE  III).   Final  glide  angle  was  minus  2 
degrees.   A  1.2  wg"  flare  was  planned. 

For  these  conditions 

From  FIG.  8 

#/,  =  /JZ3S         fa  -  32S/jpz) 

With  equation  (29)  a  stepwise  trajectory  can  then  be  con- 
structed (FIG.  13).   The  steady  glide  conditions  are  found 
to  be  as  follows: 


2? 


With  equation  (33) 


J=.    ^      J?./Z 

^_  __  jet  ft 


This  approach  was  flown  on  the  simulator  and  the 
actual  envelope  of  4  trajectories  is  also  shown  in  FIG.  1^, 
The  altitude  loss  was  900  ft.  and  T  varied  between  approxi- 
mately 13  and  1?  seconds. 

With  this  very  simple  analysis  the  ability  of  any 
particular  vehicle  to  perform  a  planned  approach  can  be 
quickly  evaluated.   Simulator  experience  has  shown  that 
an  actual  approach  can  be  very  closely  approximated  by 
this  simple  method  of  synthesis. 

The  long  period  dynamics  of  the  re-entry  glider  land- 
ing can  be  characterized  by  eight  variables.   These  are: 
L  -  Lift  V  -  Velocity 

D  -  Drag  S  -  Reference  area 

g  -  Acceleration  of  gravity  f>   -  Air  density 
tf -   Flight  path  angle        W  -  Vehicle  weight 
This  study  has  shown  that  for  typical  vehicles  the  lift- 
drag  characteristics  can  be  expressed  in  terms  of  a  para- 
bolic equation  for  the  landing  flight  regime.   If  this 
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Is  done,  the  landing  can  be  characterized  by  5  non-dimensional 
groups : 


J- 

L  I 


*'$% 


£~ 


P^e*.      SffirCc 


The  flight  path  angle  d"  varies  from  some  steady  glide 
value  given  by 

to  a  typically  small  final  glide  angle  of  about  minus  2 
degrees. 

The  normal  load  factor  "n"  is  pilot  selected  and 
essentially  constant  during  the  flare  maneuver. 

The  velocity  is  normalized  by  VR  which  contains  pure 
vehicle  parameters. 

The  maximum  lift  to  drag  ratio  E  is  a  pure  vehicle 
characteristic,  but  is  independent  of  wing  loading  allow- 
ing a  degree  of  freedom  between  the  2  vehicle  parameters 
E  and  VR. 

The  deceleration  parameter  ^  Is  a  non-dimensional 
time  in  which  the  effects  of  vehicle  parameters  on  a  de- 
sired deceleration  time  In  final  glide  are  shown. 
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PIG.  9  u  versus  T  for  final  glide,  and  FIG.  15  for 
the  final  approach  and  flare  show  how  the  5  non-dimensional 
parameters  define  generalized  trajectory  characteristics. 
From  these  generalized  characteristics  the  effect  of  the 
vehicle  parameters  E  and  Vo  can  be  shown  as  in  FIGS.  16 
and  1?.   Increasing  E  gives  lower  final  approach  glide 
angles  and  sink  rates,  less  kinetic  energy  and  altitude 
loss  in  flare,  and  longer  deceleration  times  in  final 
glide  for  a  given  initial  velocity.   Increasing  VR  results 
in  increased  final  approach  sink  rates,  greater  loss  of 
altitude  and  kinetic  energy  in  the  flare,  higher  landing 
speeds  but  a  slightly  increased  deceleration  time  in  final 
glide  for  a  given  initial  velocity. 

Since  the  vehicle  parameters  E  and  VD  are  enough  to 
separate  the  trajectory  characteristics  of  vehicles  which 
reasonably  fit  the  assumption  of  parabolic  drag  polars,  it 
seems  to  suggest  that  E  and  V^  could  be  used  as  a  basis 
for  rating  vehicles  in  the  manner  of  flying  qualities 
parameters.   If  a  pilot  did  not  object  to  the  short  period 
characteristics  of  a  re-entry  vehicle,  the  long  period 
dynamics  would  then  become  the  more  important  and  dominate 
his  opinion  of  any  vehicle.   Simulator  tests  with  Vehicles 
C  and  D  which  have  equal  E  and  VR  but  varied  wing  loading 
and  drag  polar  indicated  this  was  true.   In  this  case  no 
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significant  difference  in  vehicle  performance  was  noted. 
However,  the  experiment  was  decidedly  limited  In  this  area 
and  needs  verification.   It  seems  quite  reasonable  that 
increases  In  E  would  always  be  welcomed,  but  that  an  opti- 
mum In  VR  might  be  found  since  the  effects  of  VR  are  not 
encountered  consistently  in  all  3  phases  of  the  landing 
as  they  are  with  E. 

9.   Landing  Dispersion. 

The  average  landing  dispersion  for  the  k   vehicles 
tested  was  -300  feet  based  on  68  landings  in  which  data 
was  available.  The  maximum  deviation  from  the  mean  touch- 
down position  for  each  vehicle  was  as  follows: 

Over    Short 


VEHICLE  A 

+1000 

-900 

VEHICLE  B 
(Including  B^ 

+1000 

-1500 

VEHICLE  C 

+900 

-1000 

VEHICLE  D 

+1000 

-900 

Because  the  touchdown  and  flare  completion  speeds 
were  predictable  for  a  given  final  approach  condition,  it 
was  found  that  the  approach  glide  could  be  "aimed"  initially 
to  give  the  vehicle  a  high  probability  of  touchdown  within 
a  small  landing  area  (£1500  feet  from  a  predicted  touchdown 
point).   In  the  presence  of  turbulence  the  probability  that 
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the  airplane  might  be  "bumped"  Into  an  early  landing  In 
the  last  2  or  3  seconds  of  the  final  glide  was  found  to 
exist. 

Because  of  the  low  altitude,  and  distance  from  the 
field  (approximately  8,000-12,000  feet)  In  final  glide, 
the  pilot  may  not  have  a  clear  view  of  the  landing  area. 
For  this  reason,  It  Is  possible  that  approach  zone  range 
markers  could  be  an  Important  landing  aid  to  Improve  touch- 
down  accuracy. 

10.   Short  Period  Dynamics. 

Although  the  primary  emphasis  In  this  study  was  on 
long  period  vehicle  motions,  It  was  found  In  preliminary 
simulated  flights  that  the  basic  airplane  had  unsatisfac- 
tory short  period  dynamics  in  the  landing  flight  range. 
The  basic  re-entry  vehicle  has  of  Itself  little  or  no 
rotary  damping.   Therefore,  without  some  pitch  rate  feed- 
back, the  short  period  characteristic  motion  did  not  pro- 
vide the  necessary  tight  attitude  and  flight  path  control 
necessary  for  precise  landing.   Particularly  in  turbulent 
air,  the  unaugmented  longitudinal  dynamics  proved  unsatis- 
factory.  Moreover,  it  seems  unlikely  that  a  re-entry 
vehicle  configuration  could  perform  its  varied  mission 
satisfactorily  without  extensive  longitudinal  stability 
augmentation. 


32 


Using  a  constant  speed  approximation,  the  simplified 
short  period  motion  Is  given  by  the  characteristic  equation 

(35)      A*  +  (^f-Me-Mi)x-(M<f  L*Mi)-o 

from  REF.  6.   The  dimensional  stability  derivatives  and 
short  period  characteristics  for  Vehicles  B  and  C  are 
given  In  TABLE  II  for  speeds  of  350  and  230  fps.   Due 
to  equipment  limitations  /fy  was  Invariant  with  flight 
condition. 

The  damped  natural  frequency,  as  approximated,  showed 
fairly  good  agreement  with  test  data.   It  may  be  noted 
that  the  short  period  characteristics  for  these  vehicles 
without  the  pitch  damper  would  normally  be  rated  in  the 
"poor"  to  "bad"  region  in  the  "thumbprint"  (REF.  7).   It 
was  not  particularly  desired  to  overdamp  the  short  period 
motion  but  the  computer  limitations  dictated  fixed  values 
of  /&. 

The  damped  short  period  characteristics  were  found 
satisfactory  for  the  landing  study.   The  minor  differences 
between  the  short  period  characteristics  of  the  various 
vehicles  was  negligible  so  far  as  the  landing  problem  Is 
concerned.   The  same  was  found  with  the  difference  in 
control  power,  gust  sensitivity,  and  CL<Sf.   Future  studies 
might  consider  these  variables  in  more  detail. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  general  conclusions  are  made  for  the 
range  of  parameters  and  assumptions  of  this  report: 

1.  The  final  approach  and  landing  technique  for  the  re- 
entry may  be  divided  into  3  phases  as  follows: 

PHASE  I  -  FINAL  APPROACH.  Steady  state  glide  at  con- 
stant airspeed  and  flight  path  angle  aimed  at  a  fixed 
geographic  point  short  of  the  runway. 

PHASE  II  -  FLARE.   A  nearly  constant  load  factor  pull- 
up  during  which  the  sink  rate  is  retarded  while  the 
altitude  and  flight  path  angle  are  reduced  to  near 
level  flight  at  low  altitude. 

PHASE  III  -  FINAL  GLIDE.   Deceleration  to  touchdown 
in  essentially  one  "g"  flight. 

2.  For  the  landing  maneuver,  the  re-entry  vehicle  aero- 
dynamics may  be  approximated  by  a  parabolic  drag  polar. 

3.  For  the  parabolic  drag  assumptions,  a  re-entry  vehicle 
approach  and  landing  may  be  synthesized  based  on  5  gen- 
eralized, non-dimensional  parameters;   n,  ¥ ,   u,  E,  and  7*  . 
The  approach  is  calculated  in  reverse  from  the  touchdown 

to  the  steady  state  final  glide  matching  boundary  conditions 
between  phases. 
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**.   If  the  minimum  touchdown  speed  is  not  otherwise  limited, 
a  speed  of  approximately  .88  VR  is  the  minimum  landing 
speed  giving  acceptable  final  glide  performance. 

5.  A  minimum  deceleration  time  in  final  glide  of  10  to  15 
seconds  is  required  for  marginal  landing  performance  while 
15  to  20  seconds  yields  satisfactory  performance. 

6.  An  average  load  factor  of  1.2  to  1.3  "g"  was  pilot 
selected  as  "near  optimum"  for  the  flare  maneuver  independ- 
ent of  altitude  loss.   For  trajectory  analysis,  a  step 
function  of  load  factor  with  amplitude  of  1.2  "g"  gives 
excellent  results. 

7.  An  average  landing  dispersion  of  ^300  feet  and  a  maxi- 
mum  of  -1500  feet  from  the  planned  touchdown  position 
might  be  expected  based  on  68  landings  with  4  vehicles 
initially  positioned  on  an  optimum  final  glide. 

8.  A  maximum  sink  rate  at  touchdown  of  10  fps  is  a  prac- 
tical limitation  and  was  seldom  exceeded  for  an  optimum 
approach.   Touchdown  usually  will  occur  at  a  small  but 
finite  flight  path  angle  of  about  minus  2  degrees. 

9.  The   short  period  dynamics  of  the  basic  vehicle  are 

generally  unsatisfactory  for  precise  landing  performance. 

With  pitch  rate  damping  incorporated,  the  long  period 

dynamics  predominates  the  landing  performance  evaluation 

with  E  and  V  as  the  most  significant  vehicle  parameters. 
R 
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The  following  general  recommendations  are  made  as  a 
result  of  this  investigation: 

1.  Further  simulator  study  to  determine  the  extent  to  which 
E  and  VR  determine  a  re-entry  vehicle's  landing  performance 
as  seen  by  the  pilot;  i.e.,  make  wide  variations  of  the  wing 
loading  and  parabolic  drag  equation  parameters  while  keeping 
E  and  V„  constant. 

2.  Further  investigation  to  determine  if  T  final  glide  is 
a  function  of  E.  and  VR;  i.e.,  study  the  effect  of  decelera- 
tion, landing  speed,  and  distance  from  touchdown  at  flare 
completion. 

3.  Determine  if  an  average  flare  load  factor  of  1.2  to 
1.3  "g"  is  valid  over  a  wide  range  of  pilots  and  vehicles. 

k.      Determine  if  limits  of  E  and  V_.  are  defined  by  approach 
glide  sink  rate,  glide  angle  and  altitude  lost  in  flare. 

5.  Investigate  the  importance  of  slope  of  the  lift  curve 
on  airplane  touchdown  attitude  and  pilot  visibility. 
Investigate  the  short  period  dynamics  to  find  the  limit- 
ing pitch  damping  necessary  for  precise  longitudinal  con- 
trol. 

6.  Determine  the  effect  of  large  airspeed  and  altitude 
offsets  from  the  optimum  final  glide  conditions. 
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TABLE   I 
PROGRAMED   VEHICLE   CHARACTERISTICS 
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riTY 

UNITS 

A 

B 

Bl 

D 

5 

ft.2 

160 

136 

136 

160 

160 

M 

lbs. 

3640 

4000 

4000 

4000 

5025 

"/s 

lbs/  ft.2 

22,7 

29,4 

29.4 

25.0 

:-.4 

?> 

slug-ft? 

3698 

4520 

4520 

3698 

3698 

A 

ft. 

21.1? 

21.05 

21.05 

21.17 

21.17 

4. 

0 

+  .10 

+  .10 

-.12 

-.12 

C 

2.15 

1.725 

1.72 

2.15 

2.15 

ClS 

.814 

.573 

.573 

.814 

.814 

.0588 

.0377 

.0564 

.0575 

.0402 

(-Pot 

0 

.094 

0 

-.390 

0 

C»J 

.022? 

-.0306 

.228 

-.158 

-.0349 

Cdc< 

1.195 

1.752 

1.21 

1.49 

.391 

(-£>ot£ 

0 

2.02 

0 

+  .40 

0 

<-  rr?o 

0 

-.0072 

-.0040 

0 

0 

(-rrtc/ 

-.187 

-.0552 

-.0612 

-.187 

-.187 

Cm* 

-.268 

-.1245 

-.1146 

-.268 

-.268 

t>Y}0(& 

0 

-.0524 

0 

0 

0 

-  860 
V 

-   1260 
V 

-  1260 
V 

-   860 
V 

-   860 

NOTE:  C*   =  equivalent  rotary  damping  for  M^  10 
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TABLE  II 
DIMENSIONAL  STABILITY  DERIVATIVES 

AND 
APPROXIMATE  SHORT  PERIOD  CHARACTERISTICS 


VEHICLE  B 

VEHICLE  C 

V  fps 

V  fps 

Derivative 

230 

350 

230 

350 

.527 

.795 

.757 

1.15 

/Y& 

-10 

-10 

-10 

-10 

/f« 

0 

0 

0 

0 

Af« 

-2.2 

-5.1 

-10.8 

-25 

Frequency 

and 

Damping 
Ratio 

3 

CHe=o) 

.178 

.177 

.115 

.115 

co* 

1.^82 

2.25 

3.29 

5.0 

5 

(M&—10) 

i , .. 

2.38 

1,50 
3.6 

I.25 
^.3 

.923 
6.05 

<40 


FIGURE    I 
REENTRY    VEHICLE 


FIGURE  2 
VEHICLE  B 
LONGITUDINAL  WIND  TUNNEL  DATA 
CLEAN  CONFIGURA 
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FIGURE   3 
VEHICLE   C 
LONGITUDINAL  WIND  TUNNEL  DATA 
CLEAN   CONFIGURATION   M-.2 
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FIGURE  3 
VEHICLE  C 
LONGITUDINAL  WIND  TUNNEL  DATA 
CLEAN  CONFIGURATION  M=.2 
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APPENDIX  A 

Details  of  the  Simulation  Equipment 
A,l  Mechanization  of  Equations. 

The  longitudinal  equations  of  motion  were  simulated 
using  the  Electronics  Associates  Pace  analog  computer  and 
a  Libratrol  500  digital  computer  linked  through  the  analog 
digital  voltage  converter.  The  three  longitudinal  equations 
for  body  axes  plus  kinematic  relations  are: 

<  A-i )      ^<r=  -  £a*tr0  +jpr*g  fa  +<^«  +  4,ff 

(A-2)       ?-  -  £^<-  £V~*  **  "i/o  "^ /4b  +  <£„  * 


V-  ^>//  ^   &«s<xf  +   fa*'«*f 


(a-3)     a-  -fy y'^  Afe*,.  +  ?»,««+  ?»,,  f  +  C„/  «*j 


0 

(a-4)        A  —  iSe&t'KAx'  J-  i/s^L.x'o 

(A-5)  X=    VdO<ts^  -  i^s^^r0  A* 

(A-6)  <^-     /£-    f*-*^)  /&*«£, 

(A-7)  c^=-    £?~-^  -  4^ 
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where  the  following  approximations  have  been  made: 

<&?<L  A*  ~  J. 

The  schematic  of  the  mechanization  of  the  above  equa- 
tions is  shown  in  FIG.  Al. 

The  aerodynamic  coefficients  given  in  TABLE  I  were 
determined  by  a  best  fit  to  the  aerodynamic  characteristics 
in  the  range  of  both  trimmed  and  untrlmmed  flight  conditions 
suitable  to  the  problem. 

Dynamic  pressure  was  computed  in  the  Libratrol  ^>0Q 
using  analog  inputs  of  velocity  and  altitude  and  atmos- 
pheric density  programmed  in  the  computer.   A  standard 
exponential  atmosphere  of  the  form 

was  used  in  the  computation. 

Experimental  data  was  recorded  on  the  analog's  X-Y 
plotter  and  the  6  channel  Visicorder  906C  strip  chart 
recorder  having  a  galvanometer  natural  frequency  in  excess 
of  120  cycles  per  second.   The  X-Y  plotter  was  normally 
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employed  to  plot  trajectory  data  such  as: 
y  ys.  <^  j     /?  yx.  V;      <?/■     A  t<?  X 
The  906C  recorded  time  histories  of  the  following  variables: 


K^,  j-j   A,    X,  ^  ^</  <£ 


A. 2   Cockpit  Display. 

The  overall  view  of  the  cockpit  including'  fixed  base 
seat,  instrument  display,  CRT  runway  display,  and  side 
controller  is  shown  in  FIG.  1.   The  detailed  sketch  of  the 
instrument  panel  and  runway  display  is  shown  in  FIG.  A2. 

The  runway  display,  computed  in  the  Libratrol  500 
using  analog  inputs,  shows  a  runway  10,000  feet  by  200  feet 
in  standard  perspective.   Airplane  pitch  attitude  is  also 
an  input  to  the  CRT. 

The  y   indicator  shows  airplane  flight  path  angle  in 
degrees.   The  airspeed  indicator  was  calibrated  in  feet  per 
second.   The  attitude  indicator  read  &    in  degrees.   The  £4 
meter  was  calibrated  to  indicate  glide  path  error  in  hundreds 
of  feet.   The  direction  of  the  error  was  such  that  the  air- 
plane must  fly  toward  the  needle.   The  altimeter  was  cali- 
brated to  read  in  thousands  of  feet  to  an  altitude  of  2,000 
feet,  then  rescaled  to  hundreds  of  feet.   The  roll  attitude 
indicator  was  not  used  in  this  study.   All  indicators  were 
scaled  for  inputs  from  the  analog  computer. 
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A. 3  Computation  of  Runway  Display. 

The  runway  display  was  computed  in  the  Libratrol  500 
using  inputs  of  h  and  x  from  the  analog  computer.   The 
equations  of  the  perspective  can  be  developed  considering 
FIG.  A3. 

Computing  the  front  and  rear  of  the  runway: 


Then 


(A-8)  c£=>     -f7*tfi=*     ^(x+£) 


(A-9)  df-  J-'ioffi   -  / 


X 


Computing  the   slope   of  the  runway  edges: 

J/ 


Then 

(A-10)  ^/  =  jf   _  £ 
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Using  programmed  runway  dimensions  and  eye  to  screen  dis- 
tance, the  Llbratrol  500  computed  the  perspective  view 
with  equations  A-8,  A-9,  and  A-10.   An  Input  of  airplane 
attitude  e?  to  the  CRT  completed  the  display. 
A.*4-  Simulation  of  Atmospheric  Turbulence. 

The  atmospheric  turbulence  was  basically  a  uniform 
probability  distribution  of  angle  of  attack  with  adjustable 
RMS  amplitude  and  bias.   The  random  number  capability  of 
the  Llbratrol  500  digital  computer  was  used  to  generate  a 
square  wave  signal  F(t)  with  finite  pulses,  equal  spacing, 
random  amplitude  and  a  frequency  of  12  cycles  per  second 
as  shown  In  FIG.  A*4-(a).   The  probability  distribution  Is 
shown  In  FIG.  A4(b).   From  REF.  8  the  power  spectral  den- 
sity f (f)  is 


as  shown  In  FIG.  A4(c). 

As  shown  by  schematic  FIG.  A5(a),  the  random  signal 
F(t)  was  filtered  and  used  to  modulatea  1KC  carrier,  then 
recorded  on  a  Sony  262  C  AM  tape  recorder.   The  corner 
frequency  of  the  first  order  filter  was  1  radian  per  second. 

The  random  signal  from  the  tape  recorder  was  then 
demodulated  ,  biased  for  equal  positive  and  negative  «£  and 
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summed  with  ex:  to  give  the  desired  turbulence  level.  The 
gust  level  was  normally  adjusted  to  an  equivalent  RMS 
vertical  gust  of  6.5  fps  at  the  speed  for  minimum  drag. 
The  RC  constant  of  the  demodulator  was  .5  seconds.  The 
filtered  and  biased  turbulence  power  spectral  density  and 
probability  density  are  shown  in  FIG.  A6. 
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